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We study the throughput-delay performance tradeoff in large-scale wireless ad hoc networks. It has been shown 
that the per source-destination pair throughput can be improved from Q{l/^n logn) to 0(1) if nodes are allowed 
' to move and a 2-hop relay scheme is employed. The price paid for such an improvement on throughput is large 

^j,, delay. Indeed, the delay scaling of the 2-hop relay scheme is <d{n\ogn) under the random walk mobility model. In 

^ ■ this paper, we employ coding techniques to improve the throughput-delay trade-off for mobile wireless networks. 

, For the random walk mobility model, we improve the delay from <d{n\ogn) to 9(n) by employing Reed-Solomon 

■ codes. Our approach maintains the diversity gained by mobility while decreasing the delay. 
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I. Introduction 

In large-scale wireless networks, a fundamental question is how much information one source node can 
transmit to a destination in the presence of other simultaneous transmissions when the number of nodes 
grows large. In their seminal paper [1], Gupta and Kumar investigate the network layer throughput, or 
the capacity, of large-scale wireless networks. They study the scaling behavior of the throughput in static 
random wireless networks, where the network has n stationary nodes distributed in a unit disk. Every 
node serves simultaneously as a source for a randomly chosen destination, as a destination for another 
source, as well as a relay for all other source-destination (S-D) pairs. The authors show in [1] that when 
^ the nodes are randomly distributed, the per S-D pair throughput scales as ^Jn logn)lj and when the 
• nodes are optimally placed, the throughput scales as Q{l/^/n). This is the best achievable throughput 
performance even allowing for optimal scheduling, routing and relaying schemes. This result is somewhat 
pessimistic since as the network size grows, the throughput goes to zero. 

In [2], Grossglausser and Tse propose a 2-hop relay scheme and show that a surprisingly significant 
improvement on throughput can be achieved if nodes are allowed to move. Indeed, they show that when 
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all nodes follow a uniform, stationary and ergodic mobility model, the 2-hop relay scheme can achieve a 
constant, i.e., 6(1), per S-D throughput, which does not vanish as the network size grows. Diggavi et al. 
show in [3] that a constant per S-D throughput can be achieved even with constrained mobility models. 

The price paid for such an improvement on throughput is large delay. Characterizing the trade-off 
between the throughput and delay in mobile wireless networks has since attracted intense interest [4]- 
[11]. In [8], El Gamal et al. study the trade-off between per S-D throughput and delay in static and 
mobile wireless networks. They show that in static networks, the throughput-delay tradeoff satisfies 
D{n) = Q{nT{n)). Assuming a random walk mobility model, they show in [8] that the delay scaling of 
Grossglauser and Tse's 2-hop relay scheme is 6(nlogn). The same throughput-delay scaling behavior is 
also obtained by Lin et al. in [10] for a Brownian motion mobility model. In [5], Neely and Modiano 
study the throughput-delay trade-off in mobile wireless networks under an i.i.d. mobility model, where 
they partition the unit square into 6(n) cells and assume that after one time slot, each node moves from 
its current cell to another uniformly chosen cell. Neely and Modiano show that under this i.i.d. mobility 
model, the throughput-delay trade-off satisfies D{n)/T{n) > Q{n). Specifically, they show that the delay 
scaling of Grossglauser and Tse's 2-hop relay scheme is 6(n). A similar i.i.d. mobility model is also 
considered in [4]. Given that different mobility models and assumptions yield different throughput-delay 
trade-off behaviors, Sharma et al. propose a unified framework to capture the trade-off in a general setting 
in [9]. 

All previous work discussed thus far did not consider coding techniques. Recently, Ying et al. proposed 
joint coding-scheduling algorithms to improve the throughput-delay performance for mobile wireless 
networks using rate-less codes (e.g. Raptor codes) [12]. They study the random walk model where the unit 
torus is divided into sub-squares and at each time mobile nodes move from its current sub-square to 
one of its eight adjacent sub-squares. In this case, they show that the optimal throughput is 0{^J D{n)/n) 
when S = o(l) and D = uj{ \ \ogS\l S"^). They present a joint coding-scheduling algorithm that achieves the 
optimal throughput 0(y^D(n)/n) when S — o(l) and D is both a;(max{(log^ n)| log5'|/5'^, -y^logn}) 
and o{n/ log^ n). Note however, that with D = o{n/ log^ n), the optimal throughput is o(l/ logn), which 
is vanishing as n — > oo. Hence the algorithms proposed in [12] are restricted to cases with small delay 
and small throughput, where a constant throughput (©(!)) is not achievable. More importantly, when 
S — l/i/n, the model in [12] is identical to the one studied in [8]. However, the results of [12] do not 
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apply to this case, since the set of feasible values of D is empty. 

In this paper, we study the throughput-delay trade-off in mobile wireless networks employing Reed- 
Solomon (RS) coding. In particular, we propose a 2-Hop Relay with Reed-Solomon Coding (2HRRSC) 
scheme, and show that while maintaining a constant 0(1) throughput, 2HRRSC scheme can achieve a 
delay scaling of 9(n) under the random walk mobility model studied in [8]. The improvement of the 
delay from 6(nlogn) to 6(n) under the random walk mobility model is significant. The intuitive idea 
behind this improvement is as follows. In the 2-hop relay scheme achieving 0(1) throughput in [8], a 
particular relay node and a destination node need to meet and also be scheduled as a sender-receiver 
pair, which only occurs with a small probability. If they are not scheduled as a sender-receiver pair when 
they meet, they need to wait until they meet again and also be scheduled as a sender-receiver pair. This 
increases the packet delay. In contrast, in our scheme, the destination does not need to wait to meet a 
particular relay node, it can collect a packet from any relay nodes it encounters. This will be more clear 
as we present the detailed algorithm and analysis below. 

The proposed 2HRRSC scheme not only provides a significant improvement in the network throughput- 
delay tradeoff, but also offers potential benefits in terms of robustness and security. In networks deployed 
in challenging environments, it is possible that a relay carrying a particular packet from the source never 
reaches the destination. In the 2-hop relay scheme without coding, this would mean that the packet would 
never be delivered to the destination. In the proposed 2HRRSC scheme, however, it is not required that 
all relays carrying coded versions reach the destination. As long as a sufficient fraction (to be specified 
below) of the versions are delivered, the destination can still decode all packets. Furthermore, in the 
2HRRSC scheme, any relay carries at most a few coded versions of the packets. This means that the 
relays cannot typically decode the packets (which are intended for the destination) themselves. Hence, the 
proposed scheme is more secure than an uncoded scheme, where the relays have access to the uncoded 
packets they carry. 

The rest of this paper is organized as follows. In Section II, we outline the network and mobility models 
that we use in this paper, and define some basic concepts. In Section HI, we first briefly introduce Reed- 
Solomon (RS) codes, and then present a 2-Hop Relay with Reed-Solomon Coding (2HRRSC) scheme 
for mobile wireless networks. In Section VI we study the throughput-delay performance of the 2HRRSC 
scheme under the random walk mobility model. Finally, we conclude this paper in Section V. 
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II. Models and Definitions 

Suppose that at time 0, n nodes are unifomily distributed at random in the unit two-dimensional torus 
B. Let xi*^ denote the location of node m in ;B at time t. Each node serves as a source for one and only one 
destination, and each node serves as the destination for one and only one source. The network exhibits one 
possible source-destination pairing out of all possible pairingsjj with all pairings being equally probable. 
Since each S-D pairings is equally probable, we assume throughout the paper that the network exhibits 
a particular pairing. Let S-D pair i be the S-D pair with source node i and its destination. Assume all 
sources have saturated transmission buffers, i.e., every source node always has a packet for its destination. 

Divide the unit torus B into y^x-y/n equal-sized cellsj^and label them as {(z, j) : z, j = 0, 1, y/n—1}. 
Assume that time is slotted such that each node moves from one cell to another cell after each time slot. 
In this paper, we consider a random walk mobility model where each node performs a simple random 
walk on the ^/n x y/n grid. That is, if a node u is in cell at time t, then at time t + 1, the 

node u is in any one of the four adjacent cells of with equal probability, and the location of node 
u, Xn^^'*, is uniformly chosen at random in the new cell. By adjacent cells of we mean cells 

{{i + 1, j), (i — 1, j), + 1), j — 1)} with the addition and subtraction being modulo y/n. 

Under the network and mobility model described above, suppose a node u transmits a packet at time 
t, then a node v can receive this packet successfully if and only if for any other transmitting node w in 
the network, ||xj^^ — xl*''|| > (1 + A)||xi*'' — xi*''||, where || ■ || is the Euclidean distance and A is a 
positive number. This model is referred to as the Protocol model [1], and has been widely used in [4], 
[5], [8]-[10]. Another model for transmission is the Physical model [1]. Since these two models are 
essentially equivalent for the analysis of throughput and delay scaling performance in large-scale wireless 
networks [1], we focus on the Protocol model in this paper. 

Transmissions in our mobile wireless network are coordinated and controlled by a scheduling scheme. 
More precisely, a scheduling scheme tt is a sequence of polices {nk} which determines which nodes in 
the network transmit, and which packet is transmitted at each node, at each time slot k = 1, 2, .... For a 
given scheduling scheme, the throughput and delay are defined as follows: 

Definition 1: (Throughput) For a given scheduling scheme tt, let i?,r(^, t) be the total number of packets 

^The total number of pairings is X^fe^o(~l)'° (fc) ~ '^)'- 

'For simplicity, we ignore integer constraints throughout the paper. 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY 5 

transmitted for S-D pair i up to time t. The long-term throughput for S-D pair i is 

lim inf -HT^{i, t). 

t-^oo t 

The average throughput over all S-D pairs is 

n 

r;(n) = -Vliminf-i/,(i,i). 
n ^ — ' t^oc t 

i=l 

The throughput of tt, T^{n), is defined as the expectation over all network realizations of the average 
throughput over all S-D pairs, i.e., 

Un)=m{n)]. (1) 

Note that a network realization includes a realization of initial node positions {x|^°\ xi°^} and node 
movement trajectories {Xf\ ...,X^*^} for all t >0. 

Definition 2: (Delay) For a given scheduling scheme tt, let -D7r(i, k) be the delay of packet k (from the 
time the packet starts transmission at the source until the time the packet is decoded successfully at the 

destination) for S-D pair i. The delay for S-D pair i is 




The average delay over all S-D pairs is 




The delay of tt, DT^{n), is defined as the expectation over all network realizations of the average delay 
over all S-D pairs, i.e., 

D„{n) ^ E[D'^{n)]. (2) 

The throughput T'^{n) and delay D'^{n) are both random variables, since they depend on the initial node 
locations and random node movements. The throughput TT,{n) and delay Dt^^u) are ensemble averages. 
To study the asymptotical behavior of T-!^{n) and 0^^(11), we will let the number of nodes n go to infinity. 
We say an event holds with high probability (w.h.p.) if the event occurs with probability 1 as n goes to 
infinity. 

III. 2-Hop Relay with Reed-Solomon Coding Scheme 

In this section, we first briefly describe the Reed-Solomon (RS) coding scheme, and then propose a 
2-hop relaying scheme with RS codes that achieves 6(1) throughput and 6(n) delay under the random 
walk mobility model. 
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A. Reed-Solomon Codes 

An (n, m) linear code C over a finite field Fg is an m-dimensional subspace of the vector space of 
all n-tuples over Fq, consisting of codewords. Here q — p^is the size of the finite field Fg, where p is 
a prime number and s is an integer. An (n, m) linear code with minimum Hamming distance d is called 
an (n, m, linear code. 

A fundamental relationship for the parameters (n, m, d) over any field is the Singleton bound: d < 
n — m + 1. Any code that satisfies the Singleton bound with equality is called a Maximum Distance 
Separable (MDS) code. A useful class of MDS codes are the Reed-Solomon (RS) codes. It is known that 
for any code parameters n,m,d — n—m+1, {1 < m < n) and finite field F,, there exists a linear (n, m, d) 
(extended) RS code over ¥q, as long as n < q + 1 [13]. A natural way to define RS codes is to take 
n = q. Let the m information symbols be denoted by (t/q, Z/i, ■ ■ ■ , Vm-i) where yi ^ ¥q,Q < i < m — 1. 
Let f{x) — yo + yix + • • • + y^-ix"*"^ be the corresponding polynomial in the indeterminate x. Let 
Pi, P2, Pq be the q different elements of ¥q arranged in some order. The information polynomial f{x) 
is mapped into the g-tuple f{P2), fiPq)) over F^, where /(A) is equal to the polynomial f{x) 

evaluated at element /9j e F^: 

m—l 

/(A) = 5Z%/3f, l<i<q. (3) 

The q"^ g-tuples generated by the mapping f(x) — > {f{Pi),Pi e F^} as the polynomial f{x) ranges over 
all polynomials over ¥g of degree strictly less than m form a linear (n = q, m, d = n — m + 1) MDS 
code over ¥q. This code is called a Reed-Solomon (RS) code. 

An essential property of RS codes (or any MDS code) is that every subset of m — n — d+1 coordinates 
is an information set of C, in the sense that the codewords run through all possible m-tuples in any 
subset of m coordinates. Thus, the m information symbols can be recovered from any m — n — d + 1 
received coded symbols. This is crucial for the improvement of the throughput-delay trade-off when we 
apply this coding technique to mobile wireless networks. A series of efficient RS decoding algorithms by 
Peterson, Berlekamp-Massey, Euclid, Welch-Berlekamp and Sudan have been developed over the years. 
For details, please refer to [13]. 
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B. Coding and Relaying Algorithms 

As we explained in the network model, each node serves as a source for one and only one destination, 
and each destination has one and only one source. Consequently, there are n sessions (S-D pairs) in the 
network. In addition, each node serves as a relay for all the other n — \ S-D pairs. Due to these different 
roles, a node has different encoding/decoding and packet relaying operations. 

We treat each packet as a symbol over a finite field Fg with q = for a prime number p and an integer 
s. We require that 2' < g, where / is the length of each packet. Furthermore, it is required that n < q + 1. 
This implies that as n gets large, the size of the finite field has to grow at least with n. We assume these 
relationships hold throughout this paper. 

(i) Source Encoding: Each node u first groups every m source packets into one generation for its 
destination. In other words, for S-D pair u, each generation g = 1,2, ... consists of m consecutive 
packets as an m-tuple y"'^ = (y"'^, 1/2'^, y^^) G (Fg)*". Node u then applies the RS encoding 
algorithm to generate a "codeword" as an n-tuple (2;"'^, Zg'^? 2;^'^) ^ i'^q)"' (cf- ®). We call z]^'^ 
the kth version of generation g, k = 1,2, ...,n. The versions for one generation are stored in the 
order of version index A; in a buffer (assumed to have infinite capacity) for node u designated for its 
destination. 

(ii) Relay Storing: Each node has n — 1 buffers for the other n — 1 S-D pairs, and all buffers have infinite 
capacities. When a node w serves as a relay for S-D pair u and receives a version z^'^ for S-D pair 
u, it puts this version into the buffer queue for S-D pair u. 

(iii) Destination Decoding: As soon as the destination node of node u receives any m distinct versions for 
the same generation g, the destination node can recover the original source packets y"'^, 7/2'^, ■■■,y'^ 
by employing the RS decoding algorithm. 

In this paper, we do not consider the delay incurred by the encoding and decoding processes. 

C. 2-Hop Relay with Reed-Solomon Coding (2HRRSC) Scheme 

In our 2-hop relay scheme, we allow nodes to transmit only to nearby nodes. Precisely, a node n in a 
cell may transmit only to the other nodes in the same cell 

Definition 3: (2-Hop Relay with RS Coding (2HRRSC) Scheme) We assume that there is a hand- 
shaking protocol that allows each node to identify those nodes in its cell, the generations that those 
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nodes have already successfully decoded, and the versions which have been received for un-decoded 
generations 

When a node u is scheduled as a sender at time t (the scheduling scheme will be discussed after 
Lemma [T]), it chooses one of the following two actions: 

(i) Serves as source with probability ps'. If there is at least one node other than u in m's cell, node u 
chooses one such node f , uniformly at random, and transmits the HOL (head-of-line) version of the 
current generation for m's destination to node v. If node v is not m's destination, v performs the relay 
storing algorithm. If there are no other nodes in n's cell, nothing happens. If all n versions of the 
current generation g have been sent out, node u moves to the next generation g + 1 and applies the 
source encoding algorithm. 

(ii) Serves as relay with probability \ — ps'- If there are nodes other than u in m's cell, node u chooses 
one such node v, uniformly at random. Suppose that node v has successfully decoded generations 
1, 2, gf, then node u transmits the head-of-line (HOL) innovative version for the earliest (having the 
smallest index) undecoded generation (in this case g + l)'m m's buffer to node v. Here, an innovative 
version is a version that node v does not yet have. Node u also discards all versions for decoded 
generations 1, 2, (7, if it has any. If node u does not have an innovative version for an undecoded 
generation for v, or if there are no other nodes in m's cell, nothing happens. 

In the following section, we will show that this 2HRRSC scheme can achieve 0(1) (average) S-D 
throughput and Q{n) (average) packet delay. The improvement of the delay from Q{n \ogn) to Q{n) relies 
on the use of RS coding. To see this, note that in the 2-hop relay scheme achieving 0(1) throughput in [8], 
a particular relay node w and a destination node v need to meet and be scheduled as a sender-receiver 
pair, which only occurs with a small probability. If w and v are not scheduled as a sender-receiver pair, 
they need to wait until they meet again and be scheduled as a sender-receiver pair. This increases the 
delay for the packet carried by the node w. In contrast, in our scheme with RS coding, the destination 
does not need to wait to meet a particular relay node. It can collect a version from any relay node it 
encounters, and can perform decoding as soon as it collects any m out of the n distinct versions. 

''Protocols similar to the hand-shaking protocol discussed here are also implicitly assumed for all the networks studied in [4]-[ll]. For 
simplicity, we do not consider the overhead associated with this protocol. 
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IV. Throughput-Delay Performance of 2HRRSC 

In this section, we study the throughput-delay performance of the proposed 2HRRSC scheme. We show 
that the 2HRRSC scheme can achieve 0(1) throughput and 0(n) delay under the random walk mobility 
model. 

A. Scheduling Scheme 

We say that a node u has achieved a successful transmission if u transmits (either as a source or as a 
relay) a packet and the receiver successfully receives the packet under the Protocol model. 

Lemma 1: Let 9{n) be the fraction of nodes that can have successful transmission simultaneously at 
any given time. Then, there exists a scheduling scheme such that 



E[e{n)] = e' 



n \ n 



(4) 



where 9' = (\^/2{l + A)] + 1) ^ and A is the parameter defined in the Protocol model. 



Proof: We first show that among the n cells, there exist 9'n cells, where 9' is a positive constant 
independent of n, such that simultaneous transmissions in these cells can take place without interfering 
with each other. This is similar to the result proved in [8]. 

Consider a node m in a cell The largest possible distance between u and any other node in its 

cell is ^/2d, where d = l/\/n is the edge length of the cell. Under the Protocol Model, if all other 
transmitting nodes are away from node n's cell by a distance greater than or equal to [(1 + A)^/2~\d, then 
node n's transmission can be received successfully. Thus, as shown in Fig. [H the fraction of cells that 
can be chosen to have simultaneous transmissions is at least 

9' 4 = 1 .5. 

{\il + A)V2]d + dy (rv^(i + A)i+i)2- 

We refer to the cells which can have simultaneous transmissions as active cells. 

Given that a cell is active, a successful transmission can take place in the cell if there are at least 
two nodes in the cell. Let pi{n) be the probability that there are at least two nodes in a cell. Since the 
initial positions X^^'', X2°'*, xi°'' are uniformly distributed at random in the unit torus, the positions 
X^*\ Xg*'', X^i*-* are i.i.d. and uniform at any given time t > under the random walk mobility model. 
Thus, the number of nodes in each cell has a binomial distribution with mean value L Therefore, we have 

n-l" 



n / \ n 



(6) 
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Fig. 1. Scheduling of active cells 

Thus, the mean number of nodes that can be scheduled to have successful transmissions simultaneously 
is pi{n)6'n. The mean fraction of nodes that can be scheduled to have successful transmissions simulta- 
neously is E[9{n)\ = 6'pi{n), and we obtain (g)). As n — > cx), we have 

^0 = lim E[e{n)] = e'{l - 2e-^). (7) 

n— »oo 

□ 



Lemma [U asserts that, at any given time, a positive expected fraction E[9(n)] of the n nodes can have 
successful transmissions simultaneously. To achieve this, we need to first schedule active cells and then 
schedule transmitting nodes within the active cells. For the first part, according to the proof of Lemma [TJ 
any given cell can be scheduled to be active every jj time slots. This is the approach used in [8]. In 
fact, since the nodes are mobile, we can fix a set of cells with a regular pattern (e.g., the shaded cells 
in Fig. [T]) which are scheduled to be always active. Then, consider the following scheduling scheme for 
transmitting nodes: when there are at least two nodes in an active cell, one node, chosen uniformly at 
random, is scheduled as a sender. The sender then chooses one node among the remaining nodes in the 
active cell uniformly at random as the receiver. 

We can show that the probability that any given node has a successful transmission is also E[9{n)] 
as follows. The probability that a given node u is in an active cell is 9'. When there are k > 1 nodes 
other than u in the same active cell as u, the probability that u is chosen as a sender is rrr- Hence, the 
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probability that node u is scheduled as a sender is 

Let Y be the time between two successful transmissions for any given node. Because the underlying 
mobility process and the scheduling scheme are both Markov, we have 

= Emr 

B. Packet Delay 

We say that two nodes meet if they are in the same cell. Under the random walk mobility model, since 
each node independently follows a simple random walk, the joint position of the two nodes can be viewed 
as a difference random walk with respect to the position of one node. The inter-meeting times are just 
the inter-entering times to (0, 0) for the difference random walk on a i/n x ^/n grid. Then, for any two 
nodes u and v, the inter-meeting times are i.i.d. The following lemma provided in [8] gives the mean and 
variance of the inter-meeting time r. 

Lemma 2: [8] Under the random walk mobility model, the expected value and the variance of the 
inter-meeting time r are given by 

E[t] = n, Var(T) = e(n^ logn). (10) 

Note that Y is the time between two successful transmissions for a given node, where the role (i.e., 
source or relay) of the transmitting node is not specified. Now let Y' be the time between two successful 
transmissions for a node acting as a source. Then Y' — Yi +Y2 + ■ ■ ■ + Yh', where H' is a geometric 
random variable with parameter ps. Then 

In the 2HRRSC scheme, every source transmits one generation n times (n distinct versions). Never- 
theless, some nodes may receive multiple distinct versions and some others may not receive any version. 
The following lemma asserts that at a certain time instant after the source node sends out the first version 
of a given generation, there exists a constant fraction of nodes in the network carrying versions of the 
given generation. 
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Fig. 2. Time instants of transmission by node u, and inter-meeting times of nodes u and wi. 



Lemma 3: In the 2HRRSC scheme, let ti be the time instant when source node u transmits the first 
version of generation g. Let 



ts{n) 



n 



n 



(12) 



E[6{n)]ps PsO'piinY 

and M be the number of nodes that receive at least one version of generation g from u during [ti,t2]j 
where t2 = + ts{n) and [^1,^2] refers to the set of time slots + l,...,t2 = ti + ts{n)}. Let 

0(n) = ^. Then, for any < e < 1, we have 



Pr <p{n) < J 



1 - e 



+ 1 



e 



logn 



n 



(13) 



PsB'p'i(n) 

Proof: Let M' be the number of versions that the source u sends out during [^1,^2]- Then there are M' 
receivers receiving the versions, although the corresponding nodes are not necessarily distinct. Since in any 
time slot, the probability that node u has a successful transmission as a source is psE\6{n)\ = Ps9'pi(n), 
we have 

E[M'] = ts{n)ps9'pi{n) = n. (14) 

Suppose node u transmits a version to node wi at ti. Let Zi be the time between the ith and the (i + l)th 
meetings between u and wi, as illustrated in Fig. [21 Note that when wi and u meet again, node u may 
not be scheduled as a sender. 

Now define K = min{j G N : X]i=i > ts{n)}. In other words, K is the number of times node u 
meets wi during [ti,t2]- By definition, we have 



E 



K-l 



i=l 



< ts{n), and E 



K 



> tJn) 



Since the Zj's are i.i.d. and K is a stopping time for [Zi], by Wald's equality [14], we have 



E 



K 



.1=1 



and 



E 



K-l 



i=l 



E[Z,]E[K], 



{E[K] - 1)E[Z, 
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Since E[Zi] = E[t] = n, we have 



or 



ts(n) „r n ts(n) 

< E[K] < + 1, 
n n 



^ < E[K] < } , , + 1. 



(15) 



(16) 



Ps6'pi{n) ps6'pi{n 

Moreover, by Theorem [10] in Appendix A we havqf 

t,(n)Var(Z,) , ^^j^ejnHogn) 
Var(A) = h o(l) = — ^ h o(l) = 9(logr2). (17) 



E[Z,_ 



Let M be the number of distinct nodes {wi,W2,.-.,wm} among those M' receivers that receive a 
version from u during [ti,t2]- Let Rj be the number of times that node u meets node wj during [ti,t2]5 
where j = 1,2,..., M. By coupling methods, it is easy to show that E[Rj] < E[K]. 



Since M' = J27=i Rj, 

we have. 



E[M'] = E 



E 



■ M 

.i=i 



< E[M]E[Kl 



which yields 



n 



Psd'piin) 



+ 1 



(18) 



(19) 



Since M < n, we have E[M] = e(n). 

Now consider i.i.d. random variables Kj,j = 1, M having the same distribution as K. Note that 



A/ 



Var [Yl^^j - y^^y " ^Var(K) = e(nlogn). 
On the other hand, by the law of total variance, we have 



(20) 



M 

Var I 5^ 



E 



E 



E 



E 



M 



+ Var E 



■ M 

.j=i 



M 



Var|5^/^,|M 



+ Var(ME[ir]) 
+ Nax{M)E[Kf 
+ Var(M)G(l). 



(21) 



'Here, in the context of Theorem [Tol h = ts{n), = Var(Zi) = Var(r) = e(n^ logn) and fi = E[Zi] = E[r] = n. 
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By dlSl) and dlB, we have 

Var(M) = 0(ri log n). (22) 



Thus 



and 



E\M] 1 

E[0(n)] = > , , ^ , (23) 

Psd'pi(n) 



1 / lo£[ Ti 

Var(0(n)) = ^Var(M) = O ) . (24) 



By Chebyshev inequality, for any < e < 1, 



Pr I (i){n) < \ I < Pr (</)(n) < (1 - e)E[(i){n)\) 



^ + 1 



^ Var(0(n)) 
- e^E[(P{n)Y 
^ Var(0H) 



O 



Psd'pi{n) 

logn 



n 



This implies that (f){n) > (1 - e)Yf^ for any < e < 1 w.h.p. □ 



Lemma [3] asserts that ts(n) = „ , . time slots after the source node transmits the first version of a 

' ^ Psd'pi(n) 

given generation, with high probability, there exist (j){n)n = Q(n) nodes having at least one version of the 
given generation in the network. How these nodes and the destination node are distributed in the network 
affects how fast the destination node can collect m distinct versions to decode a given generation. Due 
to the nature of the random walk mobility model, we know that given any initial distribution, after some 
time, the distribution of each node converges to the uniform distribution. This effect can be captured by 
the mixing time. More precisely, we define the e-mixing time as follows. 

Definition 4: Let P = [Pij] be the transition matrix of a Markov chain with uniform stationary 
distribution. Let P* be the tth power of P. Define 

n 

A.(t) ^ 

Then, the e-mixing time is defined as 

Tmix{P, e) = sup inf {t : Ai{t') < e for all t' > t} . (26) 



n 



(25) 
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Lemma 4: Let P be the transition matrix of a random walk on a two-dimensional torus. Let T^ixiP, e) 
be defined as in (l26l) . then 

21ogn-l^^ ,n , /„2^ / Slogra ^2^^ 



n-l - ' ' ' - 1 _ i 

n 

Proof: Due to the symmetric nature of the random walk mobility model, the transition matrix P for each 
node is a doubly stochastic matrix. Let 1 = Ai(P) > A2(-P) > ■ ■ • > A„(P) > — 1 be the eigenvalues of P. 
Let XmaxiP) = max{A2(P), — A„(P)}. The following result gives well-known bounds on Tmix{P-, e) [15]. 



A^,,(P) log(2e 



,-1 



<T™,(P,e) < 



log n + log e ^ 



(28) 



^ — \ (p\ — ^'""v" — 1 _ x (p\ ■ 

^ '^max\^ ) ^ ^max\^ ) 

For a random walk on two-dimensional torus, Xmax{P) = ^ [16]. Substituting XmaxiP) = n ^^'^ ^ ~ n^' 
we obtain ^h. □ 



Suppose the distribution of node u among the n cells at time to is f^*") = {f[^°\ f2°\ fn'^^), then 
after t time slots, the distribution is f(*+*") = f (*o)pt^ where P is the transition matrix. The probability of 



state i is fl 



(t+to) 



ELi fk'^PL- Then we have 



(t+to) 



n 



(to) 
k 



pt 



P* - - 

n 



n 



k=l 
n 

^ E 

k=l 

= Mt) 



Hence, by the definition of TmixiP,^/'^^) and Lemma IH after cologn time slots for some constant 
Co > 3/(1 — ;^), the probability that any given node u is in any given cell is bounded between ^ 
and ^ + That is, the distribution is approximately uniform. This observation leads to the following 
lemma, for which the proof is given in Appendix B. 



Lemma 5: Suppose that ts{n) 



time slots after the source node transmits the first version of 



Psd'pi{n) 

generation g, there exist M distinct nodes having versions of generation g. Let ri{n) be the probability that 
a scheduled sender has a version of generation g after another cologn time slots where cq > 3/(1 — ^). 
Then ri{n) = 0(n) = ^ w.h.p. 

In the 2HRRSC scheme, when a relay node w is scheduled to transmit, it chooses one of the nodes 
in its cell, uniformly at random, as a destination node, and transmits a version designated for that node. 
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Similar to Lemma [H the following lemma quantifies the probability that (in any time slot) node v is 
chosen as the receiver by a scheduled sender which performs as a relay for v. 

Lemma 6: Let p2in) be the probability that (in any time slot) node v is chosen as the receiver by a 
scheduled sender which performs as a relay for v. Then 

n-2 



P2{n) = 9'{1 -ps) 



n — 1 



1 



i-i - 1-1 



n-l 



(29) 



n/ \ n ^ 

Proof: To be scheduled as a receiver, node v needs to be in an active cell, which occurs with probability 
0' . Assume there are A; > 1 nodes other than v in the active cell, then the probability that v is scheduled 
as a receiver is (l — -j^^ i. Furthermore, the probability that the scheduled sender is not the source and 
performs as a relay for node v is ^^(1 —Ps)- Therefore, we have 



n-l 



P2{n) 



fe=i 

n-l 
k=l 

e\i-Ps 



e\i~Ps) 







p 






n 






p 


.)" 




n 


n 


- 2 


n 


- 1 


n 


- 2 


n 


- 1 


n 


- 2 


n 


- 1 


n 


- 2 


n 


- 1 


n 


- 2 


n 


- 1 


n 


- 2 


n 


- 1 



1 - 



I (n-l 



k+l k\ k 



n — 2 1 — 1 
k 



n 



1 

1 - - 

n 



n—l—k 



n—l—k 



2^^ (nj 



n-l)! /I 



k=l 

n-l 



n-l 



k=l 

n 



k=2 



k)\k\ \n 



nl 



n 
k + l 

k) (n 

1 

n 



{k + l)\ \n 

k+l 



1 

k+l 



n~l—k 



1 - - 

n 



n-(k+l) 



1 



n 



n 

n—k 



n-{k+l) 



n 



n-l 



As n — > oo. 



P2 = lim P2{n) = e'{l - p,)(l - 2e-^). 



(30) 

(31) 

□ 



Based on the above results, the following lemma asserts that when m — the number of packets in one 
generation — is properly chosen, the destination node can collect m distinct versions for that generation 
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^ , 1 . Jo l . ^ 

t/n) CQlognT / t/n) 

Fig. 3. Key time instants and duration for transmission and collection analysis. 

within 6(n) time slots w.h.p. 

Lemma 7: Consider a particular S-D pair with source node u and destination node v. For any positive 
constant Ci, there exists a positive constant 5 satisfying 

^< (32) 

1 + ci + 1 

such that by choosing the number of packets in a generation to be m = 5n, node v can collect m distinct 
versions for the given generation within Cin time slots w.h.p. 

Proof: At ts{n) = ^ g/p^(„) time slots after the source node u sends out the first version for a given 
generation, there exist M = (j){n)n nodes that have at least one version of the given generation. In order 
to decode that generation, the destination node v needs to collect m versions. Let the time required for 
accomplishing this be T{n). 

As defined in Lemma [51 ti is the time instant when u sends out the first version of the given generation, 
and t2 = ti + ts{n). Consider a virtual system in which the destination node v knows ti and t2. In this 
virtual system, instead of collecting versions for the given generation as soon as possible, v waits another 
Co log n time slots after t2, where cq > 3/(1 — ^), and then begins to collect versions. Let T{n) be the 
time required for node j to collect m versions of the given generation in the virtual system. Clearly, T(n) 
stochastically upper bounds T{n), i.e., Pr(T(n) > x) > Pr(T(n) > x), for all x > 0. In the following, 
we consider this virtual system. 

Suppose wi is the first relay node that transmits successfully to v after ^3 = ^2 + cq log n, and the time 
instant when wi transmits to v h = t^^ + T[. Note that T[ is stochastically dominated by Y" , which 
is the time between successful transmissions where node v is chosen as a receiver by scheduled senders 
which perform as relays for v. By the Markov inequality and ( |29l ). for any c > 0, 

PrfF" < clogn) > 1 - 4^ = 1 - f , 

ciogn \\ogn J 

which indicates that Y" = o(logn) w.h.p. Therefore, T[ = o(logn) w.h.p. 
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Let tr{n) = cin — (cq + 1) \ogn, where ci > Since maXx.>o = tj.{n) > for any n > 0. 
The time instants ti, ^2, ^3, ti, and duration ts{n),T[, tr{n) are illustrated in Fig. [3l Let N' be the number 
of times that node v is chosen as the receiver by a relay node during [^4, ^4 + tr{n)]. Since the probability 
that (in any time slot) node v is chosen as the receiver by a scheduled sender which performs as a relay 
for V is P2{n), we have 

E[N']=t,{n)p2{n), (33) 

which implies E[N'] = &{n). 

Let K' be the number of times that Wi transmits successfully to v during [^4,^4 + tr{n)]. Then, as in 



the proof for Lemma ^ E ^f^i^ Zj 



Zi > trin), where Zj, i = 1, 2, K' , is 



< tr-{n) and E 

the time between the ith and the {i + l)th times when u'l transmits successfully to v. By Wald's equality 
and the fact that the Z,'s are i.i.d. with mean EW] = n, we have 



'-^<E[K']<'-^ + l. 
n n 

Hence E[K'] = 0(1). 

Moreover, by Theorem \W\ in Appendix A we havC 

[ciTi — (co + 1) logn]9(n^ logn) 



(34) 



+ 0(1) = e(logn). (35) 



Let be the number of distinct nodes {wi,W2, ...,wj^} among those A^' nodes that choose v as the 
receiver during [t4,t4 + U^n)]. Let R[ be the number of times that v is chosen as the target receiver by 
Wi during [^4, U + Uin)], where i = 1, 2, A^. Then E[R'^ < E[K']. Since A^' = R'i^ we have 



E[N'] = E 



E 



N 



i=l 



which yields 



E[N'[ 



< E[N]E[K'] 



E[N] > 



E[K'] 



(36) 



(37) 



Since A^ < n, we have E[N] = e{n). 

By the same argument used in (|20l) and (|2T|) . we can show 

Var(A^) = 0{nlogn) 



(38) 



"^Here, in the context of Theorem [Tol b = tr{n), = Var(Zi) = Var(r) = 6(71^ logn) and = E[Zt] = E[t] = n. 
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Each of these N distinct nodes has at least one version of the given generation for node j with probability 
r/(n). Then for any < e' < 1, if 5n < (1 - e')T]{n)E[N], 

Pr(5n < r]{n)N) > Pr((l - e)r]{n)E[N] < r]{n)N) 



> 1 



Pr((l - e')^[A^] < N) 
Var(A^) 



2 



1-0 . (39) 



n 

This implies that for any < e' < 1, if 5 < Mn)EiN] ^ ^^^^ node v can collect m = 5n versions for 
the given generation within tr{n) time slots with probability greater than or equal to 1 — O (^^)- 
By (1331), dMl) and we have 

E[N'] E\N'] tJn)p2(n) 
E[N] > -i-i > = .\\ '^^^ ! . (40) 



E[K'] - kM ^ X (kM + 1 



E[N] ^ P2C1 

cia lb — ^ k.>vj. 

Ps9'pi{n) 



Therefore, ^ > ^ as n oo. 



By Lemma Island Lemma l3l we have ri{n) > (1 — 0{^))(j){n) and (j){n) > (1 — e) ^^'^ ^'^^ 
< e < 1 with probability 1 - O (^) . Therefore, with probability 1-0 (^) 

vin)E[N] > A _ ^ Psg^Pi(n) _ tr(n)p2(n) ^ ^^^^ 



Consequently, as n ^ oo, as long as 

1 + 6*0^5 ci + 1 ' 

then with probability 1, the destination v can collect m versions for the given generation within T{n) = 
Co log n + T[ + trin) = Ciu + T{ - \ogn time slots. By the fact that T{n) <st f{n) and T[ = o(logn), 
we have T(n) < cin w.h.p. □ 



We are now ready to prove the main results of this paper. 

Theorem 8: Let Dt,^^^^ [n] be the packet delay of the 2HRRSC scheme under the random walk mobility 
model. If m = 5n for 6 satisfying (l32l) . and ci < ^7^, then DT,^^^^{n) = Q{n) w.h.p. 

Proof: Note that the packet delay is defined as the time from the instant when the packet starts 
transmission at the source until the instant when the packet is decoded successfully at the destination 
(Definition 2). 
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When a source has a generation to send, it transmits (at n time instants) n versions. This takes the 
source on average , . time slots. Hence we have D^, in) = ^l(n). 

'PsOpi{n) l'2rsc\ I \ J 

After these n versions are sent out, they are distributed in the network. All nodes other than the source 
and the destination act as relays and put the versions into their designated queues for the given S-D 
pair. Each relay node has a separate queue for that given S-D pair. Note that since the HOL versions in 
these n — 2 queues designated for the given S-D pair are not necessarily for the same generation, the 
destination does not in general collect versions according to the generation order. For instance, suppose 
the destination has successfully decoded generations 5'i,5'2, ■■■,9i-i, and starts to collect versions for gi. 
During the subsequent collecting process, the destination may collect versions not only for generation gi 
but also for generations gi+i,gi+2, ■■■■ This is because the relay nodes which have successful transmissions 
when they meet the destination and choose the destination as the target receiver may not have versions 
for generation gi. All the same, the relay nodes deliver their HOL versions, for generations gj, j > i + 1, 
to the destination. The destination can then decode generation gj,j > i, whenever it collects m versions 
for gj. 

Now consider a virtual system in which relay nodes deliver the versions only for the earliest (having 
the smallest index) undecoded generation to a destination. If a relay node does not have a version for the 
earliest undecoded generation when it meets the destination, the relay node does not transmit anything 
to the destination. Hence, in this virtual system, for the above example, the destination will not collect 
versions for generation gj, j > i + 1, until m versions for generation gi have been collected and gi has 
been decoded. Clearly, the packet delay in this virtual system is greater than or equal to the packet delay 
in the actual system. 

In the actual system, each relay node has a separate queue for the given S-D pair. In the virtual system, 
however, it is more convenient to view the n — 2 queues at the relay nodes for the given S-D pair as a 
single queue Q. In this virtual queueing system, both the arrival and departure processes are with respect 
to generations. This is illustrated in Fig. SI The packet delay in the virtual system consists of three delay 
components: the transmission delay for a generation before its arrival to queue Q (for sending all n 
versions), the queueing delay in Q, and the service time for the generation (for collecting m versions at 
the destination). 
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Fig. 4. Virtual queueing system: g,; denotes generation i. Vj denotes version j of generation i. Xi is the inter-arrival time between generation 
gi and generation Qi + l. Yj is time between the transmissions of Vj and Vj+i. Ui is inter-departure time between gi and gi+i, and Wj is 
time between the departures of Vj and Vj+i. 



For queue Q, the inter-arrival time between generations and gi+i is Xi, which is the time between 
the transmission instants of the first version of gi and the first version of gi+i. Note that Xi = J2]=i 
where Yj is the time between transmissions of two consecutive versions VJ and VJ^^ of generation gi. 



By we know E[Y'] 



PsE[e{n)] 



Thus, 



E[Xi. 



n 



(42) 



PsE[e{n)] 

Now consider a new queue Q' which has the same potential departure process as Q. The inter-arrival 
time for a generation to Q' is X" = J2j=i where Z" is the time between instants when the source 
node is in an active cell and chooses to perform as a source. In other words, in the system for Q', the 
source node sends out a version whenever it is in an active cell and chooses to perform as a source (i.e., 
it needs not be scheduled). Clearly, the queueing delay of Q' stochastically dominates the queueing delay 

of g. 

Let Z[ be the time between consecutive returns of the source node to active cells. Under the random 
walk model, since there are 9'n active cells with a regular pattern as shown in Figure [U is exactly the 
return time of a random walk on a x grid. By Lemma [21 we know E[Z[] = jr and Var(Z^) = 

6(g^log J7) = 0(1). Furthermore, Z'J = Z[ + Z'^ + h Z'jj,, where the Z^'s are i.i.d. and H' is a 

geometric random variable with parameter independent of Z^'s. Then 

1 



E[ZJ] = E[H']E[Z'^ 



(43) 
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and by the law of total variance, 

Var(Z;') = Yar{Z[ + Z'^ + --- + Z'H,) 

= E[yar{Z[ + Z'^ + --- + Z'h,\H')] + yar{E[Z[ + Z'^ + --- + Z'h,\H']) 

= E[H']Yar{Z',)+Yar{H')E[Zi]' 

= 9(1). (44) 



Therefore, 



and 



E[X':]=nE[Z';] = -^, (45) 



Var(X") = nVar(Z") = Q{n). (46) 

By Lemma 111 we know that after the destination starts to collect versions for a given generation, within 
cin time slots, w.h.p. the destination can collect enough versions to perform decoding. Consider another 
queue Q" which has an identical arrival process as Q', but with deterministic inter-service times f/f = Ciu. 
Then the queueing delay of Q" is stochastically greater than or equal to that of Q' (as can be shown by 
coupling methods). By choosing ci < y^, it is guaranteed that the generation arrival rate is strictly 
less than the generation service rate so that the queue is stable. 

By Kingman's bound [17], the queueing delay of Q" is 

which implies DT,,^^^^{n) = 0{n). Therefore we have Dt^^^^^^ti) = Q(n). □ 



C. S-D Throughput 

We now study the S-D throughput of the proposed 2HRRSC scheme. As the next theorem shows, the 
2HRRSC scheme can achieve 0(1) throughput under the random walk mobility model. The key point 
is that if each generation is transmitted Q(n) times and the generation has 0(n) packets, then each 
transmission contains = 6(1) information. Because every version is transmitted two times (2-hop 
relay), each S-D pair achieves 6(1) throughput. 
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Fig. 5. The effect of probability of being source ps. 

Theorem 9: Let TT,^^^^{n) be the S-D throughput of the 2HRRSC scheme under the random walk 
mobility model. If m = 5n with 5 satisfying (|32l) and Ci < then T^2rsci''^) = ©(1) w.h.p. 

Proof: Under the assumptions of the theorem, the proof of Theorem [8] shows that the queueing system 
corresponding to the network operating under the 2HRRSC scheme has finite 0(n) delay. For this stable 
queueing system, the throughput is the arrival rate. In our setting, the generation arrival rate is Eii^iM 
The packet arrival rate is = 5ps0'pi{n), which implies that T^2rsci''^) = Psd'pi{n)6. Since 

E[e{n)] ^ 00 as n ^ oo, we have T^^rsd^) = ©(1) w.h.p. □ 



V. Numerical Experiments 

In this section, we present simulation results on the delay performance of the proposed 2HRRSC 
scheme. In the simulation, we have n nodes uniformly distributed in the unit square, which is divided 
into ^/n X ^/n cells. The n nodes follow the random walk mobility model as described in Section II. We 
set A = V2 - 1 so that 

e' = 1 = 1 

([v^(l + A)l+l)2 9- 

This implies that each center cell of every non-overlapping 3x3 cell-grid is scheduled as an active cell. 

Figure 15] shows the impact of ps, the probability that a scheduled sender chooses to serve as a source, 
on the packet delay of the 2HRRSC scheme, for the cases of n = 225 and n = 441, and S = |. As ps 
ranges from 0.1 to 0.9, the packet delay performance follows a "U" shape, where the delays are large 
when Ps is either very small or very large. When ps is very small, encoded versions are stored in the 
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Fig. 6. The effect of parameter 5 = m/n. 

source for a long period of time, though the versions are delivered from relays to the destination quickly. 
When ps is very large, encoded versions are stored in relays' buffers for a much longer period of time, 
though the versions are sent out from the source more quickly. In the implementation of the 2HRRSC 
scheme below, we choose ps = \- 

An important parameter of the proposed 2HRRSC scheme is the encoding ratio 5 = m/n, which 
significantly affects both the throughput and delay performance. Figure |6] shows the effect of 5 on the 
packet delay of the 2HRRSC scheme with n = 225 and n = 441 nodes. From the results, we see that 
when S is relatively small (e.g., 5 < 0.5), the packet delay increases linearly with 6. On the other hand, 
when S is close to 1 (e.g., 5 > 0.8), the packet delay increases exponentially with 6. To understand this 
behavior, recall that the packet delay consists of three parts: the time required for the source to send 
out all versions, the queueing delay associated with the relay nodes (cf. proof of Theorem [8]), and the 
time required for the destination to collect enough versions to decode. When 6 is relatively small (e.g., 
S < 0.5), all three delay components increase almost linearly with respect to 6. When S is close to 1 
(e.g., S > 0.8), although the first delay component grows linearly, the third delay component follows a 
Q(nlogn) growth trend. This is reminiscent of the phenomenon in the coupon collection problem [18]. 
Also, when 5 becomes large, although the arrival rate of the queueing system of relay nodes (cf. proof 
of Theorem [8]) does not change, the service rate decreases since the destination takes a longer time to 
collect versions. This causes the second delay component, the queueing delay, to increase exponentially. 
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Fig. 7. Packet delay of 2-hop relay scheme and 2HRRSC. 

In Lemma |7l we showed that if 

^ ^ PsOq _ P2C1 
1 + 9oPs ci + 1 ' 

for 2.123 = < Ci < = 18, then the destination node can collection enough versions for a 
given generation to decode within cin time slots after the source sends out the ?7,th version of the given 
generation. Note that 0.6798 < ^ < 0.9474, = 0.5, Oq = e'{l - 2e-^) = i^f^ = 0.0294, and 
P2 = 9'{1 - ps)il - 2e-i) = i^ff^ = 0.0147. Hence 

1.4480 X 10-^ < ■ < 2.0176 x 10"^ 

1 + 6oPs ci + 1 

Therefore, if we choose 6 < 1.4480 x 10^^, we can guarantee that the packet delay has a 0(n) 
scaling behavior. Such a small 5, however, is very conservative. In order to show Lemma Ul we have 
used the conservative bound (|32l) . Choosing a small 5 deteriorates the throughput performance. In real 
implementations, 5 can be chosen to be much larger than 1.4480 x 10""^ or 2.0176 x 10^"^. For example. 
Figure |7] compares the packet delay of the 2-hop relay scheme [8] and the 2HRRSC scheme with ps = 0.5, 
6^' = I and 5 = |. From the results, we see that the packet delay of the 2HRRSC scheme follows a linear 
growth with respect to n, and also has much smaller value than that of the 2-hop relay scheme [8]. 

VI. Conclusion 

In this paper, we studied the throughput-delay trade-off in mobile wireless networks employing coding 
techniques. In particular, we proposed a 2-Hop Relay with Reed-Solomon Coding (2HRRSC) scheme, and 
showed that the proposed scheme can simultaneously achieve 6(1) throughput and 9(n) delay under the 
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random walk mobility model. The delay improvement from Q{n\ogn) to Q{n) under the random walk 
type mobility model is significant. This improvement is achieved by combining the diversity introduced 
by mobility with a reduction in delay introduced by RS coding technique. In contrast to the 2-hop relay 
scheme in [8], where the destination needs to collect any given packet from a particular relay node, in 
the 2HRRSC scheme, the destination can collect encoded packets from any of a number of relay nodes it 
encounters. In a mobile network with n nodes, the 2HRRSC scheme employs an {m,n) Reed Solomon 
code with an appropriately chosen parameter m = Q(n). The scheme guarantees that the destination node 
can decode the original m packets from the source after collecting any m encoded packets within 6(n) 
time slots. The proposed 2HRRSC scheme not only provides a signicant improvement in the network 
throughput-delay tradeoff, but also offers potential benets in terms of network robustness and security. 

Appendix A 

Theorem 10 (Theorem 5 in [19]): Let Xi,X2, ... be i.i.d. with mean /i > and finite positive variance 
cr^. Let 5*0 = 0, Sn = XliLi -^i define r(6) = mi{n : Sn > b}. Then 

Var(r(6)) = ^ + — + 0(1), (48) 

where K is a constant. 

Appendix B 

Proof of Lemma\5\- By Lemma HI we know after cq logn time slots with cq > 3/(1 — ;i-), the probability 
that any given node u is in any given cell is bounded between Ps = — ^ and P4 = ^ + 

Senders are scheduled uniformly at random in each active cell. Consider an active cell («, j). r]{n) is 
the probability that the sender scheduled in cell has a version of generation g. Let ki be the number 
of nodes within cell that have versions of generation g, and let k2 be the number of nodes in 
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that do not have versions of generation g. Then 

M n-M 

Vin) = ^Pr(A;i = i)Pr(A;2=j) 

I -r J 



i=l j=0 
M n-M 

^ EE 

i=l j=0 

M n-M 



EE: 



i=l j=0 

n min{fc,M} 
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M\ fn-M 



= E E i f rf)^'(--'"- 
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E^4(l - ^3)" V ^ ( ^ - 1 j ( ^ - ^ J + ^ ^^^^ - T ^ ( ^ - 1 
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i=0 
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.M fn 



fe=i ^ ' 

0W[(i+P4-P3r-(i-P3y 



0(n) 



1 + 



0(n) 1 + 



n 



where we used Vandermonde's identities Eto C^^"') (^1"",) = (feli) foj" ^ < and E.=o ' C^^"') (^1""^) = 

(^ij) for A; > M+ 1. These identities can be justified by expanding (1 + + a;)"~^ = (1 + 
By the same approach, we have 

r]{n) > 0(n)[(l+p3-P4r-(l-P4n 



1 - 



1 1 
_ ( 1 ^ 



0(n) 1-0 
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Therefore, ■q{n) — (l){n) w.h.p. □ 
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